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The motion of a particle is studied in a random space-time metric, 
using a non-relativistic approximation. The randomness induces a 
' diffusion of the particle in coordinate space. Hence it is shown that 

the evolution of the probability density of the particle's positions is 
given by Schrodinger equation. 
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^ 1 Introduction 



Einstein never accepted quantum mechanics as a fundamental theory of na- 
ture. Furthermore, he beUeved that a fundamental theory could not be found 
5^ I starting from the current form of quantum mechanics (e.g. by adding hid- 

den variables) but within a completely different framework, probably that 
of general relativity. In his own words: "I do not believe that quantum 
mechanics will be the starting point in the search for this basis just as one 
cannot arrive at the foundations of mechanics from thermodynamics or sta- 
tistical mechanics" jT] . The present paper attempts to explore the possible 
derivation of a fundamental theory of motion in agreement with Einstein's 
expectations. Here I shall restrict the study to motion with small velocity in 
order to obtain a non-relativistic approximation. As we shall see, although 
the starting point is different, the theory has some similarity with the de 
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Broglie-Bohm hidden- variables theory (or Bohmian mechanics) j2] . In partic- 
ular, both theories assume the existence of trajectories for the particles. The 
difference is that Bohmian mechanics rests upon the hydrodynamic interpre- 
tation of Schrodinger equation, in which particle trajectories never cross each 
other. In contrast, in this paper I consider that the particle's motion consists 
of a random motion superimposed to a smooth one. Thus we might speak 
about an aerodynamic interpretation, which makes the approach similar to 
stochastic mechanics 

General relativity starts from the assumption that space-time may be 
curved. The curvature can be derived from the space-time metric, once a co- 
ordinate system is defined. The specific assumption in this paper is that the 
metric is random. There are several reasons for this hypothesis. Firstly, noise 
is quite natural in our very complex universe, therefore to assume the exis- 
tence of randomness is more plausible than to assume its absence. Secondly 
quantum theory contains a random ingredient in the form of vacuum fluctu- 
ations. At a difference with the standard assumption of nineteenth-century 
physics that randomness is always associated to finite (nonzero) temperature, 
quantum physics of the twentieth century contains "zeropoint fields", that 
is some randomness even at zero Kelvin. Thus I propose that the motion of 
bodies should be always studied in a random space-time metric. The ran- 
domness would be specified by defining the probability distribution in the set 
of possible metrics, but in the present paper I shall not state that distribution 
and use only some assumptions about it. 

2 Motion in a random metric 

I consider a particle which is placed in xiat time tiand in X2 at time t2- It 
moves under the action of a potential U(x) in a space-time metric 

ds\ = gooc^dt^ + 2gojcdtdx^ + gjkdx^dx^^ j, k = 1,2, 3, (1) 

where goo^goi and gij are functions of (A; ^,t), c being the velocity of light. As 
appropriate for the non-relativistic approximation to be introduced later, we 
distinguish the time, t, from the space coordinates, {x^,x'^,x^} , these defin- 
ing the position vector, x. The randomness is taken into account assuming 
that there is a probability density, -P (A) , A G A, in the set A of space-time 
metrics. (For clarity in the physical arguments to be introduced later, I use 
the notation P (A) dX for the probability distribution, without any claim of 
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mathematical rigour.) Some ambiguity appears due to the freedom existing 
in general relativity for the choice of the coordinate system, but this ambi- 
guity should not produce any confusion in what follows. In any case I shall 
assume that the choice of coordinates for every metric is made so that the 
potential U{x.) has the same functional form in all of them. 

For given A, that is a fixed metric, the motion of the particle would 
correspond to the minimum of the action 

A{X) = r {-Mcdsx - U (x) dt} = ^ L (A; x, x,t) dt = min, (2) 



L (A; X, x,t) = -Mc^ gooc^ + 2gojCXJ + gijX^x^ ~U (x.) , (3) 

where are the velocities (time derivatives of the position coordinates of 
the particle at a given time) and f/(x) is the potential. The generalization to 
include forces not derived from a potential is straightforward, but it will not 
be considered in the present paper. An alternative procedure to Q would 
be to start from the differential equation of motion for the particle in a given 
metric (i. e. specified by a value of A) but a variational principle is a more 
convenient starting point for my purposes. From now on I will speak about 
the set of Lagrange functions Q, rather than the set of metrics, but use the 
same label. A, for both sets, which should not introduce any confusion. 

In this paper I will consider only motions such that x^ << c, and conse- 
quently I shall use a non-relativistic approximation. Thus I may replace (0) 
by an expansion to second order in the velocities x^ , which gives 



L(A;x,x,t) = -M 



^ ^ / _92k 9oj9ok \ ^j^k 



(4) 

For later convenience I define the "mean Lagrange function Lq" , which might 
be obtained by an average over A of (0]), for fixed x, x and t. It is 

Lo (x, x,t) = J dXP (A) L (A; x, x,t) = -mc^ + ^mx^ - (x) . (5) 

The simplicity of eq.© is a consequence of the invariance properties which I 
assume for the probability distribution of metrics, that is invariance to trans- 
lations, rotations, time-translations and Galilean transformations. These in- 
variance properties allow for the renormalized mass m to be different from 
the bare mass M. The Lagrange function (0)) may be written as a sum 

L (A; X, x,t) = Lo (x, x,t) + Li (A; x, x,t) , (6) 
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such that all randomness goes in the second term. 

The probability distribution in the set of Lagrange functions induces a 
probability distribution in the set of possible paths of the particle, which 
may be formalized as follows. We consider the set, A^, of paths going from 
(xi,ti) to (x2,t2), each path specified by the equation of motion x-'(z/;t), 
u & N, t & [^1,^2] ; where we assume that x-' as a function of t, for fixed z/, 
possesses continuous second derivatives. Our problem is to find the proba- 
bility distribution in corrresponding to that in A. Now the condition (0) 
defines a function, u = g{X), associating a path to every Lagrange function, 
which may be got from the variational problem 

1= f P{X) S{u-g{X)) dg{X) /*' L (A; x (z/, t) , x (z/, t) ,t) rft = min, (7) 

J A Jti 

where -P (A) > (see • In fact, any function g (X) which does not 
associate to every A the path giving the minimum value to the action A (A) 

will lead to a value of J (|7j) which is greater (strictly, not smaller) than 
the one given by the function leading to the minimum value of I. 

In order to devise a practical method to solve the variational problem 
((Tj) I suppose that the time integral may be approximated by a sum over a 
finite number of times, that is e, 2e, 3e, ... with e small enough. For this we 
require the position, y, and the velocity, w, corresponding to every path at 
every one of the discrete times. Thus we write, instead of ((Tj) , the variational 
condition 

J = [ P (X) 6 {u - g (X)) dg (X) K = min, (8) 

J A 

K = J'' dt J d^y J d^w 6{y - x (z/, t))6 (w - x (z/, t)) L (A; y, w,t) , 

where 6 () is the three-dimensional Dirac's delta and the unknown function 
to be found is g (A). Now for any g (A) (not necessarily the one providing the 
minimum value of J ) we may define a probability distribution, /g(y,w,t), 
in the phase space of positions and velocities via the integral 

/,(y,w,t)= / P(A) H^-gW) dg{X) 5(y-x(z/,t))5(w-x(z/,t)). 

J A 

(9) 

Now our problem is to get, from the variational condition (jH)), another one 
involving the phase-space distribution /g(y,w, t), rather than the function 
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g, if this is possible. The use of the phase-space distribution leads to a fluid- 
dynamical picture of the motion where the actual particle is replaced by a 
statistical ensemble of particles whose phase-space density is /(y,w,t). 

A case where the variational problem (jHJ may be easily written in terms 
of a phase-space distribution corresponds to L not depending on A, that is 
when the Lagrange function is not random. I shall solve this rather trivial 
case in the hope that it may provide a clue for a more general method to be 
developed in the next section. In this case inserting Q in (jSJ leads to the 
variational condition 

dt J rf^x J d^w /(x, w,t) L(x, w,t) = min . (10) 
Putting the Lagrange function (0) in this variational problem gives 



J^' dt J d^:si J d^w /(x, w,t) 



-mw^ — U fxl 



(11) 



where we ignore the first term of (0) which, being a constant, is irrelevant in 
the Lagrange function. Without any loss of generality we may write 

/(x, w,t) = p (x,t) /;,(x, w, t), 

with the conditions 

J h{x, w, t)d^w =1, J wh{x, w, t)d^w = v (x,t) , (12) 

so that p {^,t) is the density of particles at (x,/!:) , and v (x,t) is the mean 
velocity of those particles present at that space- ime point. Now I shall make 
the minimization of in two steps. In the first one I will search for the 
function /i(x, w,t) making the action a minimum with given p (x,t) and 
V (x,t) . In the second step I will find the equations of motion of these two 
functions. 

The minimum of pip with fixed p and v requires, for every (x,t) , 

J h(x,w,t) d^w =min, (13) 

with the constraints (fT^ , which leads to 

/i(x, w,t) = 5 (w — v(x, t)) , 



5 



S {) being the 3-dimensional Dirac's delta. This means that all particles 
present at (x,t) possess the same velocity, that is the phase-space probability 
becomes 

/(x,w,t) =p(x,t)5(w-v(x,t)), (14) 

which is a phase-space dependence typical of hydrodynamics. Thus density 
and velocity are related by the continuity equation 



dp 

— + Vj =0, J =pv, 



(15) 



j being the current density. As a conclusion of the first step, the form ()14|) 
allows to state the variational condition in the form 



dt J rf^x p(x,t) 



-mv (x,t)^ — U (x) 



mm, 



(16) 



with the constraint (fTK|) . 

Now we derive the equations of motion of p and v from (|TT)|) , which is 
straightforward. We introduce the condition ()15|) in the variational problem 
by means of the Lagrange multiplier S'(x, t) and get 
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dt j d^l^ \p 



-mv^ — U (x) 



S 



dp 
di 



+ V (pv) 



0, 



(17) 



where the variation of the action must be zero for any independent variations 
of p, V and S ( the variation of S just reproduces (US))). The function S must 
fulfil the condition 

S(xi,ti) = S(x2,t2), (18) 

relating the initial and the final times and positions. From the variation of 
V in (jl7|) we obtain, after an integration by parts. 



m 



(19) 



(The integrated term is zero because p vanishes at infinity). If (|T9|l is inserted 
in (|T7jl we get, after two appropriate integrations by parts. 



5 ^ ' (it J d^yip 



,,, ^ 9S 



0. 



(20) 



6 



The time integration by parts requires that the integral 



^(x,t)p(x,t)rf3x 



has the same value at times ti and t2, which holds true in view of (|18|) . In 
the variational problem (j^UI) the variation of p leads to the Hamilton- Jacobi 
equation whilst the variation of S gives again the continuity eq. fll5|) , taking 
(fT^ into account, that is 



These results reproduce standard equations of motion in analytical mechanics 
as it should. 

3 Derivation of Schrodinger equation 

When Li ^ 0, that is the metric is random, it is not obvious that the problem 
of the motion may be solved in terms of the two functions p (x,t) and v (x,t) . 
Indeed, the condition ()14|) will not be fulfilled in general. Nevertheless, I shall 
assume that such a solution of the problem is still possible provided that we 
use a functional of p (x,t) and v (x,t) , rather than a function as in . 
In any case the randomness of the metric will produce randomness in the 
velocities of the particles arriving in position x at time t. We shall take 
into account that randomness modifying both the continuity eq. ()15|) and the 
variational condition (|T6|) . If w is the velocity of a particle placed in x at 
time t, the position of the particle at time t+ At will be, neglecting terms 
of order At^ and higher. 



Ax being a random displacement. I propose to treat the displacement 
Ax as deriving from a white noise independent of time, position, and initial 
velocity, thus leading to the following probabihty distribution for the (vector) 
displacement Ax 




(21) 



y = X + wAt + Ax, 



(22) 




(23) 
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This assumption rests upon the hypothesis that the change of position in- 
duced by the space-time dependence of the coefficients, g^^, of the metric is 
more rapid than the change induced by the external forces (deriving from 
the potential U{x.) ), which is plausible for a non-relativistic theory as ours. 
Actually, if (j^^ is correct, there is a finite (small) probability that the ve- 
locity of light is surpassed, but this fact is not a real problem within the 
non-relativistic approximation. 

From it is straightforward to get, from the density, p (x, t) , of parti- 
cles at time t the density at time t + At. We get 

p (y, t + At) = Jd'^l d^'w /(x, w ,t)Q(y - x - wAt), (24) 

where /(x,w,t ) is the phase-space probability distribution at time t. To first 
order in At the integration is straightforward if we use the expansion 

/(x,w,t)~/(y,w,t)+^(xj -?/j)— + - ~ yi)(^k - Vk)- 



(25) 

Inserting ()25|) in ()24j) and performing the integrals in x and w we get, after 
some algebra, 

1^ + Vj =0, j =pv - D Vp, (26) 
where we have taken into account 

J /(y, w ,t)d^w =p(x,t ), 1 w /(y, w ,t)d^w =pix,t )v(y,t), (27) 

(see (jl2p .) Thus we arrive at the following 

Proposition 1 The density, p(x,t), and the mean velocity, v(x,t), of the 
particles in the statistical ensemble representing the actual particle, fulfil de 
continuity eq.( 



We see that now there is a "diffusion current" (last term) in addition 
to the "hydrodynamical current" (compare with (|T5|) ). In order to get the 
substitute for |TI)|) we introduce the random velocity u = w — v(x,t), so that 



w^)^^ = v(x,t)^ + (u^)^^ + 2v(x,t)- (u) 
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where ()^^ means average over those particles present in x at time t. In this 
equation the random- velocity average, (u)^^ , is the diffusion velocity which, 
from , should be 

The difficult problem is to calculate the random square mean velocity, (u^)^ ^, 
which is not possible without a detailed knowledge of the probability distri- 
bution of metrics. Thus I will make the most simple assumption, namely 
that it is a constant independent of x and t. The lack of a clear foundation 
for this hypothesis is certainly a weak point of the present derivation, which 
will be studied more carefully elsewhere. In any case, once the assumption 
is accepted it is straightforward to arrive at the following 

Proposition 2 The mean kinetic energy of the particles ( of the statistical 
ensemble) present at (x, t) is given by 

r(x, t) = (v2 - 2Z}p- V ■ Vp) + To, 

where To is a constant. 

The assumption made in , that the change of position induced by 
the space-time dependence of the metric coefficients is more rapid than the 
change induced by the external potential, implies that the constant Tq is 
rather large so that the kinetic energy is always positive. As a conclusion the 
action (|T?)|l should be replaced by 

dt J d^^L = min, L = ^m (pv^ - 2Dv ■ Vp) - p?7(x) (28) 

where we have ignored the constant Tq which is irrelevant in the variational 
problem. Now we follow the same steps leading from (jlSp and (fTBj) to (j^UI) . 
We introduce the constraint (|26|) in the variational problem (j^Hj) using a 
Lagrange parameter S, which leads to 



Sj^'^dt J rf^x Il + S 



^ + V{pv- DVp) 



=0. (29) 



Hence the variation of v gives, after an integration by parts, the following 
relation 

,^0^ + ^. (30) 

p m 
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When this is inserted in (j26|) we get the continuity equation in the form (|2ip. 
On the other hand when (|3Up is inserted in (j29j) we obtain 



*^ ' 3 



ti 



5 dt d'^ { - p 




2m ^ ' 2 \ p 




0, 



(31) 

where we have ignored two terms whose sum equals the divergence of the 
vector field DS'Vp. We must assume that this vector field vanishes at infinity 
whence the integral of its divergence is zero. 

It is remarkable that, although we have started from two expressions, 
()26|) and ()28|) , neither of which is invariant under the reversal of time, the 
Lagrange density in (jHT|) is time-reversal invariant in the sense that the 
operation (t —t) is equivalent to just changing the sign of the auxiliary 
function S, which does not change the physics. In (j3T| it is easy to see that 
the variation of p leads to 

±.^VSf,Ui.).§^2,nD^^^0, (32) 

where we have assumed that the space-time integral of d{Sp) jdt vanishes, 
a hypothesis already made in the previous section (see below eq.(j2ni))- As is 
well known the continuity eq. ipTj) and the dynamical eq.flH^ may be obtained 
by separating the real and imaginary parts in the Schrodinger equation 

= -—v^m + U^. (33) 

dt 2m ^ ' 

provided that we identify 

h = 2mD, ^ = v^exp ("y^ . (34) 



4 Discussion 



The derivation of Schrodinger equation given in this paper might be consid- 
ered a step in the direction of Einstein's expectations, as commented in the 
introduction. Alternatively it may be seen as just a new derivation from for- 
mal assumptions (the two propositions of the previous section) devoid of any 
deep physical meaning. I cannot argue too strongly in favour of the former 
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possibility, and this is why I include the word "heuristic" in the title of the 
paper. 

In our derivation the "wave-function" \l/ is just a mathematical func- 
tion, with range in the complex numbers, whose modulus gives information 
about the probability density of the particle position and whose phase is 
related to the particle velocity via eq. (j3(Jj) . In comparison with the stan- 
dard Hamilton- Jacobi equation, see (^01) , eq. (jH^ contains a term deriving 
from the fact that the space-time metric differs from the standard one of 
Minkowski space. Thus the metric plays the role of the "guiding wave" in 
the de Broglie-Bohm theory It is not contrary to the intuition that this 
"wave" may modify the scattering cross sections, making them different from 
those derived from classical mechanics. Also it is easy to understand the ex- 
istence of stationary states of electrons in atoms, as a balance between the 
attraction by the nucleus and the diffusion caused by the random metric. 
The picture is here similar to that provided by stochastic electrodynamics |3] 
with the gravitational field (the non-Minkowskian metric) substituted for the 
electromagnetic radiation. More difficult is to understand intuitively how the 
non-Minkowskian metric may give rise to the observed sharp spectral lines of 
atoms or to the interference fringes in two-slit experiments. If the latter two 
effects derive from the metric, the mechanism is not clear from the derivation 
here presented. 

A question which arises is whether the assumption of a random metric is 
really necessary for the two hypotheses introduced in the previous section. 
Indeed, derivations of Schrodinger equation from the hypothesis of a random 
motion, not involving the space-time metric, have been recurrently proposed 
during the last 50 years or more[3],j3]. In my view there are three reasons 
why the randomness of the metric is relevant in the derivation given in the 
present paper. The first one is the fundamental character of the space-time 
background, as emphasized in the introduction. The second reason is the 
existence of a diffusion in coordinate space rather than in velocity space, as 
would be more natural in any approach starting from Newtonian mechan- 
ics. Also, if the origin of the random motion is gravitational, dissipative 
effects may be negligible. The third reason is that metric randomness makes 
plausible a formulation starting from a variational condition, due to the fun- 
damental role which geodesies play in general relativity. Nevertheless the 
fact that Planck's h is a. universal constant barely follows from the existence 
of random metrics. Indeed within general relativity (or Newtonian gravity) 
the motions (of small particles) are independent of the mass, which suggests 



11 



that the diffusion constant D (see (j34j) ) , rather than should be universal. 
I have no clear response to this objection. I may mention only the fact that, 
when we consider the diffusion of a system consisting of n particles of equal 
mass, m, it may be shown that the diffusion constant of the center of mass 
of the system is D/n if the particles diffuse independently. In this case the 
product, mD, of the mass times the diffusion parameter is the same for the 
individual particles and for the center of mass of the system, independently 
of n. 

An interesting question is whether the derived equations ((EH) and (jS21)) 
are equivalent to Schrodinger equation. The answer is in the negative. The 
reason is that the derivation of the previous section implies that the function 
S (x, t) must be single- valued whilst in Schrodinger theory the wavefunc- 
tion \& (x, t) is single-valued, which just requires that S'(x, t) changes by an 
integer multiple of 271% along any closed line. Actually both assumptions 
are equivalent, by continuity, if the region of definition of S (or is sim- 
ply connected. However it is currently assumed that there are experimental 
situations where this is not the case. In particular the popular two-slit exper- 
iments, which have already been performed with electrons, neutrons, atoms 
and even molecules, are currently analyzed assuming the existence of a region 
forbidden to the particles, namely the one occupied by the screen with the 
slits, so that the region allowed to the particles is not simply connected. How- 
ever it might be possible to study these experiments replacing the screen by 
a high, but finite, repulsive potential so that the region is simply connected. 
Another case where the region of interest is not simply connected corresponds 
to some excited states of atoms. In this case the wavefunction may become 
singular at the origin (the position of the nucleus). For instance the wave- 
function of the states with quantum numbers Z = m 7^ contains the factor 
exp{il(j)) and the function S changes by 2iTlh in a rotation by 27r. A related 
problem is the existence of nodal surfaces in some solutions of Schrodinger 
equation, that is surfaces where p = but v 7^ 0. These solutions are un- 
physical in our approach. In summary every physical solution of ()21|) and 
(j221) is a solution of Schrodinger equation , but there are solutions of the 
latter which are not solutions of the former. This seems to imply that the 
formalism here developed cannot agree with the empirical evidence. Never- 
theless I question the current wisdom that all solutions of the Schrodinger 
equation are really necessary for the interpretation of the experiments, but 
this point will not be analyzed further here. 

In any case the formulation here presented allows for an interpretation of 
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Schrodinger equation in terms of trajectories, which may be useful in some 
apphcations. In this sense the formulation is similar in spirit to Bohmian 
mechanics|2]. However our trajectories present a random element and so 
the picture achieved is actually more similar to stochastic mechanics. But 
stochastic mechanics is currently understood as fully equivalent to Schrodinger 
theory, which leads to counterintuitive behaviour like the existence of nodal 
surfaces mentioned above. For this reason it is usually considered as just 
a formal approach to quantum mechanics rather than a different physical 
theory 13.^. 

The formalism here presented might be extended to many-particle sys- 
tems by replacing the three dimensional space by the 3A^ dimensional con- 
figuration space of N particles. However all particles will move in the same 
space-time metric, which would induce correlations in the motion, in addition 
to those derived from the possible inter-particle forces. Thus the generaliza- 
tion is not trivial. Incidentally, I guess that the said correlations might be 
related to Bose statistics, but this point will not be discussed further here. 

In summary, the formalism here developed allows an interpretation of 
Schrodinger equation in terms of particle trajectories, which may have some 
interest. However there are great difficulties to take it as a physical theory 
underlying quantum mechanics. 
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